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1 Introduction 

One of the most difficult problems in the theory of the Navier-Stokes equa- 
tions are related to a stationary two dimensional flow in an exterior domain, 
namely to the problem 

v ■ Vv - uAv + = F in ft, (1.1) 

V-v = in ft, (1.2) 

B(v) =0 on dtt, (1.3) 

lim v(x) = v^, (1.4) 

\x\— *oa 

where the sought solution (v,p) is a velocity vector field and the correspond- 
ing pressure, v is a viscous positive constant coefficient, F - an exterior force 
acting on the fluid, - a prescribed constant vector field and B{v) stands 
for a boundary conditions, e.g. Dirichlet boundary conditions, as v — v* on 
dfl. In our case the system will be supplemented with the slip boundary 
conditions, namely: 

n ■ T(v,p) ■ t + f(v • t) = b on dCl, (1.5) 
v ■ n = on dCl, (1.6) 



where / is a nonnegative friction coefficient, T(v,p) is the Cauchy stress 
tensor, i.e. T(v,p) = uD(v) + pi, where D(u) = {v it j + Vj^}? - =1 is the 
symmetric part of the gradient Vt> , and I is the identity matrix. Moreover 
n, f are respectively the normal and tangential vector to boundary dfl of an 
exterior domain Q, i.e. Q = R 2 \ B, for a bounded simply-connected domain 
BCK 2 . 

The slip boundary conditions govern the motion of particles at the bound- 
ary - relation (11.51) is just Newton's second law. From the physical point of 
view this constraint is more general than the Dirichlet boundary data, since 
for / — > oo and b = one can obtain relation vsn = 0. The case where / = 
is important for applications, since then the fluid reacts with surface dfl as 
the perfect gas ([3], [T5]). 

In many modern applications, like the model of motion of blood, polymers 
and liquid metals, this type of boundary conditions is widely used ([6], [T2]). 
Our considerations in an exterior domain are also important for example in 
the field of aerodynamics, where problems with flow past an abstacle is of 
high interest. 

There are many questions related to this problem, namely: existence 
of solutions, uniqueness and asymptotic behaviour. In this paper we are 
concerned with the first issue for arbitrary data. 

We are concerned with weak solutions to the problem (jl.ip - (jl.4D thus it 
is natural to require finite Dirichlet integral 

Vv\ 2 dx < oo. (1.7) 

Since the power 2 coincides with the dimension of the domain we are not able 
to use standard embedding theorems for v to get some information about the 
velocity at infinity and assure that (11.41) holds in any sense. Many mathe- 
maticians brought their attention to this problem. Some partial results were 
obtained by Gilbarg and Weinberger ([9], [10]), like the case with B(v) - 
Dirichlet boundary condition with = on dQ. Then the following asser- 
tions hold: a) every solution to (11.11) - (jl.3p that satisfies (jl.7p is necessarily 
bounded; b) for every solution to (II. ip - (11.31) that satisfies (jl.7p there exists 
Voq such that 

i>2tt 

lim / \v(\x\,8) - VoofdO = 0. (1.8) 
M^ 00 Jo 

In 1988 Amick pQ published a paper, where he proved that if the body B is 
symmetric around the direction of i^, and boundary data is symmetric 
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with respect to the direction of v^, then there exists a symmetric solution 
v,p such that 

lim v(x) = Voo, uniformly. (1-9) 

\x\^oo 

These results however give no information about the relation between 
and Voo (see [8] for more detailed information about this problem). 

On the other hand Finn and Smith [5] and Galdi [7] showed that for 
small values of Reynolds number A and when ^ there exists at least 
one solution to the system (ll.ip - (ll.4l) in a proper space. This has been done 
by applying contraction mapping technique and proper //-estimates for the 
Oseen system in exterior domain. In our paper we would like to point out 
that this technique should also work for our problem considered with slip 
boundary conditions. In [16] we give proper L p -estimates for the exterior 
Oseen system with slip boundary conditions. 

The main purpose of this paper is to show that the system (|l.ip -( fOj) 
together with (ll.5l) - (ll.6p admits at least one weak solution for arbitrary data. 
This result is gatherned in the following theorem: 

Theorem 1.1 Let v > 0, f > 0, F e (Vii 2 (fi))* and b G H-V 2 (dn). Then 
for a properly constructed vector field vq there exists at least one weak solution 
(in the sense of Definition \2. 1\) 

v = (v oo + v ) + V ± <p (1.10) 

to the system U. ij) -( f7T"7| ) with boundary conditions B(v) as in U.5\) - l[J~b]) . 

for which the following inequalities holds: 

h\\H$M<C = C(v,f,b,F) and \\Vv\\ L , {n) < C(u, f, b, F). (1.11) 

Note: We used here a simplified notation VHq(Q) which stands for: 

Vif 2 (fi) = {V<^ : <p e Hl(tt)} C Hl{n). (1.12) 

See also (2Tty . 



While dealing with a problem fll.ll) - fll.4p one faces a problem of a kernel 
of the rot-div operator. In bounded simply connected domains this kernel is 
trivial and full information about velocity v can be retrieved from its vorticity 
a. However in case of unbounded non-simp ly connected domain some more 
precise way of finding solution should be followed. In Section 11.11 we show 
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such results for bounded non-simply connected domains, namely we show 
that the kernel of this operator, when a solution satisfies the slip boundary 
conditions, is trivial. 

Our construction of v is done with the Galerkin method, what might be 
useful for the process of finding the solution numerically It has one more 
attribute - it can be easily modified to approximate a solution with solutions 
considered in a sequence of bounded domains, for which we are able to show 
that their kernel part of the mentioned operator is trivial. 

That is why we may conclude that the cons tructed vector field is a 
solution with a trivial kernel part of the rot-div operator. However we may 
not exclude the case, that there exists a solution to our problem with a 
nontrivial kernel part. 

The result from Section 11.11 can be applied to improve results from [H] , 
namely the results from this paper can be extended (with proper assumptions 
on the boundary data) to the case of non simply-connected bounded domain. 

One of the classical approaches to show existence of solutions is to use 
the Hopf inequality (see [TTj ) 

< C\\u\\ H i {n) for all^Gtfo 1 ^), (1.13) 

where 5(x) = dist (x, dfl) is the distance from the boundary, to get a'priori 
estimates on a solution u. However, in case of slip boundary conditions we 
are not able to get full information about the velocity on the boundary dfl 
(in particular u ^ Hq(Q)), as it is the case for Dirichlet constraint. Thus 
different methods are needed to get a'priori estimates. 

We introduce a reformulation of our system in terms of the vorticity 
a = rot v of the fluid. This increase the order of equations and transforms 
slip boundary conditions on v into Dirichlet condition on a. This way will 
be cleared up in the Section [21 

Our paper is organized as follows: in Section 11.11 we give some results 
about the rot-div operator and a reasoning how to improve results from [H] . 
In the next section we give a reformulation of our problem in terms of rotation 
of the fluid, introduce a function space Hq(Q) and give a definition of a weak 
solutions to our problem. The next section is the main part of our paper, 
where we obtain a'priori estimate for our sought vector field. This estimate 
allows us to show existence of solutions via the Galerkin method (see: Section 
EJ). 
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1.1 The kernel function ip. 

Note: We want to emphasize, that the following considerations concern a 
bounded domain Q, while the domain for the main problem fll.ll) - fll.4p is an 
exterior one. 

We would like to discuss some details related to the kernel function ip 
for the rot-div operator. The kernel function ip appears when one wants to 
get back the information about velocity vector field v from its rotation a. 
Namely, the following system is considered: 

rot v = in Q, (1.14) 
V-v = infi, (1.15) 
v ■ n = on dCl. (1-16) 

From fll.15p - fll.16p and from the Poincare Lemma one can take v in the form: 

v = V ± tp, (1.17) 

for some function ip, which we call the kernel function. From (11.161) on can 
see that 

ip = const on dQ. (1.18) 

Since rot V ± if) = Aip we can figure out that ip = (up to an additive con- 
stant) in bounded simply-connected domains. This was the case in |T4]. In 
this section we would like to show that the assumption on simply-connectedness 
of the domain in [13] can be ommited. 

Before that we want to point out that the family of solutions to the system 
fll.14p - fll.16p . in the case of a bounded domain with the first homotopy group 
IIi(O) = Z, is one dimensional. This is a simple observation, since all these 
solutions can be presented in the form 

V = C^, (1.19) 

where is a suitable constant and \1> is the unique solution to the system: 

= in Q, (1.20) 
^ = onTi, (1.21) 
* = 1 on T 2 , (1.22) 

where and T 2 are two disjoint connected parts of dQ. 
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In [13] it is shown that also for an exterior domain a family of solutions to 
a similar problem is also one dimensional. In that case, however, one needs 
to use additional property for a solution v, namely v — > as \x\ — > oo. In 
the same paper a constant is explicity given for a linear modified Oseen 
problem. 

What is interesting is that the slip boundary conditions allow us to cal- 
culate that the kernel part of the solution is equal to zero. Indeed, in our 
problem the slip boundary conditions under rot operator transforms into 

Totv = (2x-f/v)v-T ondQ, (1.23) 

thus from (11.141) we see that 

(2* = 0. (1-24) 

Since it is impossible that (2% — f/v) = on all dQ (because of positivity of 
/ and v) it has to be 

at least on an open subset of dQ. This is however impossible since \I/ is a 
harmonic function which takes its maximum and minimum on the boundary 
dQ, and from the strong Maximum Principle 

M > < L26 > 

at points, where \I> equals its maximum, which contradicts (jl.25p . The case 
where \I> = const in the whole Q stays in contradiction with boundary con- 
ditions fQ31) - fTCT|) . 

This conclusion might be used to extend results from [14J to the case 
where the domain is non simply-connected. In this paper authors are consid- 
ering simply-connected domains because in this case the kernel of the rot-div 
operator is obviously trivial and one may recover full information about the 
velocity of the fluid from its rotation. 

2 Reformulation 

In this section we give a reformulation of our problem in terms of the vorticity 
a of the velocity vector field (see [T9]). 
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Taking rotation of (II. ip we get: 

- uAa + v ■ Va = rot F, (2.1) 

where a = rot v = v 2> i — t>i,2- The slip boundary conditions give us full 
information about the vorticity of the fluid on the boundary, namely from 
fll.7l) - fll.8p we get a condition on a: 

a = (2 X - f/v)v -r + b on dn, (2.2) 
v -n = 0, (2.3) 
v — > Voo as |x| — * +oo, (2.4) 

where x is the curvature of the boundary. The exact procedure was consid- 
ered in [19|. 

Next, we take over the information at infinity from v. Let us introduce 
an extension vector field v = v + satisfying the following conditions: 

V • v = in fi, (2.5) 
v = on dtt, (2.6) 
vo — > Voo as |x| — > oo. (2.7) 

The vector field vq will be defined later. The construction will fulfil require- 
ments from Lemma [2.21 in order to get a'priori estimates on ip. 
Having v we rewrite v as 

v = v + u, (2.8) 

where for u we have the following constraints: 

V -u = in A, (2.9) 
u ■ n = on90, (2.10) 
m — ► as \x\ — > oo. (2-11) 

Having fl2.9p - fl2.10l) we use the Poincare lemma to present u in the following 
form: 

u = W ± ip } (2.12) 

where V ± (p = (—^,2? Since u ■ n = |p = on <9f2 we may take <£> = on 
<9Q. 
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For futher calculations let us notice that 

rot = Aip. (2.13) 

We may now derive from (jl.ip - fll.4D the system of equations for u. Re- 
calling that 

a = rot v = A<p + rot vq (2-14) 

we write: 

-i/Aq + (v + VV) ■ Va = rot F in fi, 

a — Atp = rot Vq in Q, 

tp = on dfi, (2.15) 

a - (2 X - //^)V> ■ f = (2 X - //i/)t5o • ^ + & onffl, 
V" 1 ^ — >■ as |x| — >■ oo. 

Since we want Dirichlet integral for v to be finite, i.e. 

|2 



Vv\ 2 dx < oo, (2.16) 

we establish a suitable space for the solution (p. Let us introduce the following 
Banach space: 

if 2 (f2) = C °°(fi) l|V2 ' lli2(n) (2.17) 



equipped with a norm 



^ |VV| 2 ^ • (2.18) 

We now define a weak solution to the problem fll.lD - fll.4p . 

Definition 2.1 We say that v is a weak solution to problem to- 
gether with fli. 7| )-( TJ7^) iff there exists p> G Hq(Q) such that v = v + V -1 ^ 
and the following identity holds for every ip G Cq°(Q): 

/ (V"V + v ) • Vip{Atp + rot v )dx + v I (Atp + rot v )Aip dx+ 
Jn Jn 

-v f (2x-/H(VV + tF )-Aa= fF-V^dx. (2.19) 

To prove existence of solutions in the sense of Definition 12.11 we first need to 
show a'priori estimates on <p, what will be the case in Section 12.11 Before 
we get into this, we would like to investigate some properties of the kernel 
function of the rot-div operator, to explain the properties of the constructed 
solution. 
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2.1 A'priori estimate. 

To show existence of a solution to our problem we use the standard Galerkin 
method. This is a construction of approximate solutions which converge in 
some sense (strong enough to pass to the limit in the equation) to the limit 
vector function. The construction of this sequence requires usage of a fixed 
point theorem. Proper converging of this sequence (or a subsequence) can 
be assured by showing uniform boundedness (in a proper function space) of 
all its elements. In both of these steps a great help is a'priori estimate of 
solutions to our equation. That is why we now focus on obtaining it. 

We follow the standard approach to get a'priori estimate, i.e. we multiply 
(I2.15li ) by if and integrate over Q. Recall that 

a = Aip + rot v . (2.20) 

The first term —u j n Aaip dx gives us: 

— v I Aaip dx = —v \ (A<p + rot v ) Aip dx+ 
Jn Jn 

vf (2 X -fM((v + V ± ip)-T + b)^;da. (2.21) 
Jdn an 

For the second term it is easily seen that 



v + V ± (f) ■ Vapdx = (v + VV) • Vyxt dx (2.22) 
n Jn 

since (p = on dQ and V • (v~o + V -1 ^) = in fi. Finally, since V ± (p ■ Vy? = 
we may write 

/ (v + VV) • Vaipdx = / v ■ V<padx. (2.23) 
Jn Jn 

This term causes difficulties in getting a'priori estimates for the solution. 
The reason why is that it is of order 2 with respect to <p (see 12.201) . just like 
J n \Aip\ 2 dx, but without any information about its sign. Thus we need to 
prove an inequality in the form: 



Vq ■ VipAifdx 



<7l|A^||£ 2(n) , (2.24) 
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for some constant 7, which should be small enough (in our case v/2). This 
is done by a proper construction of the vector field v~ . 

Since the construction of the vector field vo is done in a neighbourhood 
of dQ we introduce compactly supported Vq such that 

v = v + v oo , (2.25) 

with proper constraints for Vq, which will be precised in the following lemma: 



Lemma 2.2 For every e > there exists compactly supported Vq, which 
satisfies the following conditions: 



V • v = in Q, 

v = —Voc on dQ, 



(2.26) 



and the following inequality holds 



(vo + Voo) ■ VpAip dx 



< e||A^||| 2(n) (2.27) 



for every ip G Hq{VL). 

Proof of the Lemma 12.21 First we transform the term 

/ vo-V^A^dx, (2.28) 
Jn 

using integration by parts, to get a term without Voq in it. This is because 
the term 

/ Voo ■ VtpAtpdx = I (p^Aipdx (2.29) 
Jn Jn 

could cause a great difficulty in estimate, since a'priori we do not know 

whether or not (p t i G L 2 (Q). 

Remark: in the following calculations we ommit boundary integrals, 

since vq = on dQ, i.e. Vq = and v^ + v^ = 0. We calculate f!2.28j) as 

follows: 



d ■ V(fA(pdx = I (vq + Voo) ■ V(fA(pdx (2.30) 
v^(p^A(p + / (^ 2) + VoojtppAtp dx (2.31) 



=: h + h (2.32) 
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h = (^o 2) + ^00)^,2^,11^ + / ( u o 2) + v oo)p,2¥,22dx =: hi + I22 (2.33) 
Jn Jn 

Recalling that vq = on dQ we calculate: 

hi = - (v^ + Woo), if,2f,idx - / {vf 1 + ^00)^,21^,1^ (2.34) 
Jn 

= - / v^}(p, 2 (p,idx + o / ^ol^A- ( 2 - 35 ) 

z Jn 

Similarly 7 2 2: 

^22 = ^ ^(^o 2) + ^oo)(v^ 2 2 ),2^ = ~J v$ip 2 2 dx. (2.36) 
For h we use the fact that vjp = on <9fi: 

Ji = / ufV.iV.nda: + / ufV, lV&dx =: I u + I 12 . (2.37) 
Jn Jn 

Now we calculate as follows: 

In = \ I vPirthdx = ~f v^dx (2.38) 



and for 1\ 2 : 



I 12 = - I v$(p tl (p t2 dx- f v^ip^ip^dx = - I v i ^ip il ip t2 dx+^ I v^ip 2 2 dx. 
Jn Jn Jn * Jn 

(2.39) 

Finally, summing up all above calculations we get: 

I = ~ J v^tp^dx - I v$(p,i(p, 2 dx + - J v£}(p%dx 
z Jn Jn / Jn 

~ I v $. ( P,i ( P,2dx+ - / vjQp^dx- 
Jn / Jn 



I (!) , (2)\J 



Jn 


(2.40) 


/ ^(fi^dx 

Jn 


(2.41) 


v$)dx 


(2.42) 
(2.43) 



n 
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Now, since V • v = we may write: 



I = ~ <P*iVo}dx- / V 2 2 vfldx- I <p,i<p,2(v$+vj$)dx (2.44) 
Jn Jn Jn 

= - Vip- Vv ■ Vipdx (2.45) 
Jn 

In this form we see, that there is no term with i^, but its structure does 
not allow us to go into more subtle analysis of its behavior near the boundary 
of the domain. That is why we transform it into more appriopriate term. This 
is not straightforward since simple calculation by parts would lead us to the 
point we have started with. This is because there is still information about 
Voq in this term - it occurs in vq on the boundary. Thus an auxilary vector 
field is needed to take away i>oo from the boundary. We proceed as follows: 

I — — I Vy?-Vw -Vy? = - / Vcp-V(v +V e )-Vcp+ / V^-W e -V^, (2.46) 
Jn Jn Jn 

where V e is contructed as follows: let us introduce a vector field V in (t 1 ,t 2 ) 
coordinates (see Appendix): 

Wl ,y ): .(W + «)^ (247) 

for which the following conditions are valid: 

Vipih, 0)) = , V(p( tl , 0) = Q , (2.48) 

where ( = C(^) is a constant from the construction of the mapping p(t\, ^)- 
The similar conditions are fulfilled by the vector field 

V £ (p{ti,t2)) - { Q for(e<t 2 <(. (2 ' 49) 

From the construction of V e it is easily seen that 

||W e || La(n ) <C(n,V)4p- (2.50) 
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Indeed, from the definition (I2.49P we calculate: 



I |VK| 2 < C(p) [ [ 6 \V t {V(p{t x ,t 2 /z)))\ 2 dt 2 dh (2.51) 
Jn Jo Jo 

<C(p) f [ iVtVipit^h/e^l^dhdh (2.52) 
Jo Jo e 

= C(p) [ L f | W| 2 -<^ < 9M).\\ VV f (2.53) 
Jo Jo e e 

We may now estimate the integral /: 

I = - V<p- V(w + K) • V(^c/x + / V(p-VV e - Vipdx =: h + I 2 . (2.54) 
Jn Jn 

Since v o + V € = on dQ we may integrate I\ by parts and get: 

h= f Vcp(v + V e )Acpdx + l- [ ((v + V e ) -V)(\V<p\ 2 )dx = I n + h 2 (2-55) 
Jn 2 J n 

Now, from the Stokes theorem and since V ■ v q = in Q 

ii2 = -~ / (V-K)(|Vyf)cfe. (2.56) 
2 Jn 

Gathering all these calculations we get: 

I=f Vpv Apdx+ [ V(pV e A(pdx - - f (V • V e )(\ V^\ 2 )dx+ 
Jn Jn ^ J n 

+ / • VK ' Vy?cfe =: Jj + J 2 + J 3 + J 4 . (2.57) 
Jn 

Before we estimate these integrals let us introduce the following notation: 

fl e = {x e O : dist(ac, <9fi) < eC}. (2.58) 
Integrals J3 and J4 are similar and we estimate them first. Since supp V e C 

J3 + J4<c([ \VV e \ 2 ) / ([ |VyfV (2.59) 
\Jn e J \Jn € J 
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We use the following interpolation inequality: 

c\\Vv\\ LH n e) < IMI^nJI VVllJSJn.)- ^ 60 ) 
Since <p = on <9f2 and 9? G H 2 (Q e ) we use embedding theorem 

# 2 (a) C a < 1 (2.61) 

to conclude that 

Mmn t) < C(a)||VVlU 2( n e) • e 1 ^. (2.62) 

Indeed: 

|2 ^ I \,„U j- M2| 



IMl2 W </ / M*i,*2)IV;p|<fMt2 (2-63) 
JO Jo 

<C(fl) f J e t^WVMl^dtidt, (2.64) 

C(^)l|Vvii 2 (o^ 1+2a - (2-65) 



< 

Inserting this inequality to (I2.60P we get: 

l|Vv>|l*<n.) < C7(n)||VVHL(n.) • (2.66) 
Now from (12301) . (12391) and (12361) we conclude: 

J 3 + ^4 < C7(n, VOc^ll VVllL(n)- (2-67) 

Since < a < 1 (in particular a can be taken a > 1/2) we may choose e 
small enough to get 

^3 + ^4 < |||VV||| 2 (Q)- (2-68) 
The estimate of the integral J2 is similar: 

J 2 = f VipV e A<p< ( [ \V e \Ux) ' ( [ \V<p\ A dx) ' ( [ |VV| 2> ) 7 
in \Jn e J \Jn t J \Jn J 

(2.69) 

and since 

(^J \V t \ A ' < (J \Vf d^ ' =C(V)< 00 (2.70) 
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we may combine (I2.66P with (I2.7(jp and (12.691) . to get, for e small enough, 
desired estimate: 

(2.71) 



h < -||vV" 2 



L 2 (n) 



For integral J\ we refer the Reader to Lemma 2.1. from [14], which states 
that for every e > the vector field vq can be constructed in such a way that 
the following inequality holds: 



(f ■ V(p) 2 dx 



< e\\A^\\ 2 for all <p G H*(Q). 



(2.72) 



Thus integral J\ can be estimated using Schwarz inequality and above lemma: 



1 

4 



' VtpvoAipdx 
u 

This completes the proof of Lemma 12. 2[ □ 



< TllvVll^nj 



(2.73) 



Above inequalities allow us to get a'priori estimates. Namely ip fulfills 
the following identity 

-v\ \A<p\ 2 dx + v f ((2 X -f/u)((v + V X ip)-r)+b)^da+ 
Jn Jan on 

+ Vq • Vv9(rot Vq + A(p) dx = v I (rot Vq)Alp dx + / F ■ V -1 ^ dx. 
Jn Jn Jq 

(2.74) 

We see, that there is still a terms of order two with respect to <p, for which 
one cannot verify its sign, namely 

We deal with this problem using the following identity (see [13] ) for v = vq+u: 



a 2 dx 



Ian 



a(v ■ f) da 



(v) dx+ ((f • r) 2 f - b(v ■ f)) da, 



(2.76) 

which comes from the well known identity for v G H l (Q) with V • v — in 
Q (see [21]): 



a 2 dx 



(v)\ 2 dx + 
15 



'an 



(2.77) 



Since 



-v J \Ap\ 2 dx+is j ((2x-fM({vo+V ± (p)-T)+b) f||J da- J A^rot v dx = 

— 1/ / a 2 c/x+z/ / a(u-r) <icr+ / Ay?rot ■uo cfe+ / (rot vq) 2 dx—v / a(uo-f) da, 
Jn Jan Jn Jn Jan 

(2.78) 

we may derive from (12.741) using (I2.76P the following identity: 

■ / D 2 (V^ + v Q ) dx+ 
Jn 



— v 



+ / At^rot v dx + / (rot v ) 2 dx 
Jn Jn 

-v [ ((V ± ¥ > + v )-r) 2 f-b((V ± ^ + v )-r)da- (2.79) 
Jan 



—v \ a(vo ■ f) da + / t> • V(f(rot Vq + A<p) dx = / F • V <£> 
Jan 7n Jn 

To get a'priori estimate from (I2.79P we need to use the Korn's inequality (see 
Lemma 14.11 in Appendix) : 

t\2 / V7-L, „\ J™ \ I IV72,„|2 



V ± ip)dx > K / \V 2 <p\ 2 dx, (2.80) 
n Jn 

where K is a constant dependent only on Q, which allows us to get J n \V 2 (p\ 2 dx 
in the estimate, namely: 

-v I B 2 (V ± (p)dx-u [ ((VV + w ) ■ r) 2 f < -Kv [ \V 2 p\ 2 dx. (2.81) 
Jn Jon Jn 

Here we also used the fact that / > 0. Combining (12.791) with (I2.8ip we are 
able to get an estimate of the form: 



/ \V 2 V \ 2 dx < C(DATA) ( f |VV| 2 
Jn \Jn 



1/2 1 



Vq ■ VipAifdx 



n 



(2.82) 



This is because \J n vq ■ 'VipAipdx\ is the only term of order 2 with respect to ip, 

1 /2 

and all other terms are of order 1 can be estimated by C(DATA) ( f„ | V 2 y?| 2 ) 
(see the Remark below). 
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We now use Lemma 12.21 with e = Kv/2 to estimate remaining term of 
order 2 and get the following inequality: 

HVVHi^nj < C(DATA)||V^|| L2( n), (2.83) 

where in C(DATA) one includes all constants dependent on Q, F, u, etc. 
This inequality gives us of course a'priori estimate on || V 2 y?IU 2 (f2) : 

||VVlU^(fi) < C(DATA). (2.84) 

Remark: It is not hard to estimate terms in (12.801) . where u (i.e. V ± (p) 
appears in a form different than V 2 (p - one must recall that ip = on dQ, 
which gives us, together with ip e H 2 oc (Q), the fact, that all local estimates 
(the only needed) can be obtained using || \7 2 ^ || ^ 2 (f7) . 



3 Existence. 

In this section we use standard Galerkin method to prove the existence of a 
solution to 

- uAa + v ■ Va = rot F, (3.1) 

vn = 0, (3.2) 

rot v = a, in Q (3.3) 

V • v = 0, in tt (3.4) 

a = (2x - f/v)v -r + b ondtt (3.5) 

in the sense of distributions, i.e. in the sense of Definition 12.11 This means 
that we prove the existence of a function ip e Hq(Q), which satisfies (12.191) 
for all ip G C°°(Q) with compact support in Q. As was mentioned before, we 
seek for a solution in the form 

v = v + V V (3-6) 

Since Hq(Q) is Hilbertian and separable we take a base of compactly 
supported functions {lUi} 



oo . 



H 2 (n) = { Wl ,w 2 ,...r lH ^K (3.7) 
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Next we introduce a finite dimensional subspace V N (Q) C Hq(Q): 

V N (Q) = span{w 1 ,...,w N }. (3.8) 

We additionally assume that (iVi,Wj) V N = where (-,-)v N is the inner 
product in V N , which comes from the inner product in Hq(Q). We search 
for an approximation (p N of the function cp in the form: 

N 

V N (x) = J2cf Wj eV N . (3.9) 

j=0 

To find coefficients we solve the following system: 

/ (V ± ip N + v ) ■ Vwi(Aip N + rot v )dx + v (Aip N + rot v )Awi dx+ 
Jn Jn 

-v f ((2 X -fM(y ± ( P N + v )'T + b)^da= (F-V^Widx (3.10) 
Jon on J n 

ioTi = l,...,N. 

Let us introduce a mapping P : V N — > V N as follows: 

P( V N ) = y] ( f (V-V + v ) ■ Vwi(A<p N + rot v )dx 
i=i ^ 

+z/ / (A<p N + rot w ) Aw* rfx 
in 

-i// ((2r/MvV + «o)'?+t)^ 

Jen ^ 
— y F • V ± Wi dx^j ■ Wi. 

From the definition of the mapping P we easily calculate that: 

(P(0, = / ( V± ^ + • V^(A^ + rot v )dx 

Jn 

+u I (Aip N + rot v ) Aip N dx 
Jn 

-uj ((2* - /MVV + v ) ■ r + b) ^-da 
F ■ V L <p N dx 

n 



To get proper estimate for the term (P(<p N ), tp N )v N first we must use the 
following identity (see flTTSJ) or (12777ft ): 



v 



a 2 dx- / a(v ■ r)da = / W(v) dx + ((« • r) 2 / - 6(v • f))dcr. 
n ian in Jan 

(3.11) 

Hence, since V -1 ^ • V<p N = 0: 

(P{ip N ),<p N ) vN = [ v ■ Vip N {Aip N + rot v )dx - v [ amtv 
Jn Jn 

[ T> 2 (v)dx + [ ((v ■ r) 2 / - b(v ■ f))j da - f F -V L ^> N dx (3.12) 
sin Jan J Jn 

In this form it is easily seen, that one can obtain the following estimate 

(P(<p N ),<p N ) vN > for > k (3.13) 

for some constant k > 0. One must just repeat the reasoning from a'priori 
estimates. 

Such condition gives us (see [15]) existence of ip^ such that H^Hyjv < k 
and moreover P(ip N ) = 0, which solves our problem for coefficients . Thus 
we get sequence of approximating solutions <p N G V N such that H^Hyjv < k 
for some constant k > independent of N. 

Passing to the limit. Since we have uniform bound on Hv^llij^), i.e. 
1 1 #2 m) < k ( we can take a subsequence which is weakly convergent to 
some limit. However, for the sake of passing to the limit in nonlinear terms 
of (12.191) we must use diagonal technique. 

Let us denote Qr — ft fl Br. In bounded domain Qr we have tp N G 
H 2 (Qr), since tp N , Vip N G L 2 QC (Q). Moreover, since (p N = on dfl we have 

(3.14) 

hence we may choose a subsequence (f N ' , which we futher denote for simplic- 
ity as <p N , which is convergent on Qr to (p in the following sence: 

V 2 <p N VV weakly L 2 on Q R , (3.15) 

V<f N -> Vip strongly L 2 on Q R , (3.16) 

(f N — > (p strongly L 2 on Qr, (3-17) 

(p N -»• 93 strongly L 2 on <9fi. (3.18) 
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We repeat this treatment for R — > oo to get subsequence tp , which is con- 
vergent in above sense on all bounded domains. Since a test function in 12.191 
has compact support we may pass to the limit. Thus <p is a solution in the 
sense of Definition 12.11 

4 Appendix 

(ti, ^-coordinates 

We introduce (t\, ^-coordinates as follows. Let s : [0, L] — > R 2 be a 
normal parameterization of boundary dQ, i.e. 

s([0,L]) = dfi, s(0) = s(L) = x e dQ, and \s'(t)\ = 1 (4.1) 

for a fixed point xq and L - the length of dQ. Next we introduce the following 
mapping p : [0, L] x [0,Q "~ > ^ 2 sucn that 

p(ti,t 2 ) = s(ti)+t 2 n(s(t 1 )) > (4.2) 

where n is the inner normal vector to boundary dQ. If £ is small enough 
(comparing to curvature \ of boundary dQ), then the map is one-to-one and 
p G C 1 . Moreover 

dist(p(t u t 2 ),dtt) = t 2 . (4.3) 
Using the definition we compute the gradient of map p as follows 

P, h = (1 + t 2 xMs(ti)), p A = n(s(ti)). (4.4) 

Then we see that 

P,i-Lp,2 and (Vp) _1 = (— — — f,n) T . (4.5) 

By coordinates {ti,t 2 ) we denote coordinates obtained using mapping p. 
Korn's inequality 
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Lemma 4.1 For an exterior domain Q C M 2 , which is not spherically sym- 
metric, there exists a constant K > 0, dependent on the domain Q such that 
the following inequality 

/ D 2 (w) dx>K \Vu\ 2 dx (4.6) 
Jn Jq 

holds for every u G C^°" V " i2(n) satisfying 

V ■ u = in fi, u ■ n = on dVl. (4.7) 

See |2T] for the proof of this lemma. 
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